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Abstract 

Using linear-response density-functional theory, magnetic excitations in the striped phase of 
CaFe2As2 are studied as a function of local moment amplitude. We find a new kind of excitation: 
sharp resonances of Stoner-like (itinerant) excitations at energies comparable to the Neel temper- 
ature, originating largely from a narrow band of Fe d states near the Fermi level, and coexist with 
more conventional (localized) spin waves. Both kinds of excitations can show multiple branches, 
highlighting the inadequacy of a description based on a localized spin model. 
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Magnetic interactions are likely to play a key role in mediating superconductivity in 
the recently discovered family of iron pnictides; yet their character is not yet well under- 
stood. In particular, whether the system is best described in terms of large, local magnetic 
moments centered at each Fe site, in which case elementary excitations are collective spin 
waves (SWs) called magnons, or is itinerant (elementary excitations characterized by sin- 
gle particle electron-hole transitions) is a subject of great debate. This classification also 
depends on the energy scale of interest. The most relevant energy scale in CaFe2As2 and 
other pnictides ranges to about twice the Neel temperature, 2X^^40 meV. Unfortunately, 
neutron scattering experiments have focused on the character of excitations in the 150-200 
meV range l|, l2|, much larger than energy that stabilizes observed magnetism or supercon- 
ductivity. Experiments in Refs. , while very similar, take completely different points 
of view concerning the magnetic excitations they observe. There is a similar dichotomy in 
theoretical analyses of ~t,e .nteraofonsB. Q. Model deset.pfo,. usually postulate a 
local- moments picture [5[. Most ab initio studies start from the local spin-density approxi- 
mation (LSDA) to density functional theory. While the LSDA traditionally favors itinerant 
magnetism (weak on-site Coulomb correlations), practitioners strongly disagree about the 
character of pnictides; indeed the same results are used as a proof of both localized and 
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itinerant descriptions [3|, 1^ . 

The dynamic magnetic susceptibility (DMS), is the central quantity that uniquely char- 
acterizes magnetic excitations. It can elucidate the origins of magnetic interactions and 
distinguish between localized and itinerant character. However, the dynamical linear re- 
sponse is very difficult to carry out computationally; studies to date been hmited to a few 
simple systems. Here we adapt an all-electron linear response technique developed recently 
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to calculate the transverse DMS, x(q, w). Besides SW excitations seen in neutron 
scattering, we find low energy particle-hole excitations at low q, and also at high q. In 
stark contrast to conventional particle-hole excitations in the Stoner continuum, they can 
be sharply peaked in energy (resonances) and can be measured. 

The full transverse DMS x(i"7r',q, w) is a function of coordinates r and r' (confined to 
the unit cell). It is obtained from the non-interacting susceptibility xo(r, r', q, w) via the 



standard relation 
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X = Xo[l-Xo/]"' (1) 
I is the exchange-correlation kernel. When computed within the time-dependent LSDA 
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FIG. 1: (a) N{E), in units of eV~^ per cell containing one Fe atom. Data are shown for M=OAfiB 
and l.lfiB- Both kinds of transitions are reflected in peaks in xo{q=0,lj), shown in panel (c) for 
M=0.4, 0.8, and l.l/is- (b) N(E) at the Fermi level Ep as a function of moment M. The blue 
vertical bar denotes the experimental moment, and also approximately demarcates the transition 
from itinerant to localized behavior. N(Ef) in the nonmagnetic case is shown as a dotted line, (c) 
bare susceptibility xo(9=0,u;) in the same units. The text discusses the significance of the arrows 
in panels (a) and (c). 
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(TDLDA) / is local: /=J(r)5(r — r') [8|. Xo can be obtained from the band structure 
using the all-electron methodology we developed . To reduce the computational cost 
we calculate Imxo and obtain Rexo from the Kramers Kronig transformation j^. Imxo 
originates from spin-flip transitions between occupied states at k and unoccupied states at 
k + q: it is a /c-resolved joint density of states D decorated by products P of four wave 
functions Q 



D(k,q,a;) = f{el){l - fiet^^^Mu - 6^^, + el) 
Imx-^/ci.ci3kP(r,r',k,q)xD(k,q,.) 



(2) 
(3) 



Even with the Kramers- Kronig transformation, Eq. poses a huge computational burden 
for the fine frequency and k resolution required here. We make a simplification, map- 
ping Xo onto the local magnetization density which is assumed to rotate rigidly. The full 
Xo(r, r', q, 1^) simplifies to the discrete matrix xoij^,'^',^,^) associated with pairs of mag- 
netic sites (R, R') in the unit cell; and /(r) simplifies to a diagonal matrix Jrr. In Ref. [3] 



we show that we need not compute / exphcitly but can determine it from a sum rule. / 
can be identified with the Stoner parameter in models. We have found that for Fe and Ni 
these approximations yield results in rather good agreement with the full TDLDA results, 
and expect similar agreement here. We essentially follow the procedure described in detail 
in Ref. [3], and obtain x(q, cu) as a 4 x 4 matrix corresponding to the four Fe sites in the unit 
cell. To make connection with neutron experiments, spectra are obtained from the matrix 
element XIrr' < ^^'^'^Ixi^i R-'? q, Ci;)|e*'' "^' >. For brevity we omit indices RR' henceforth. 

As we will see, the character of Xo{^., ^) can largely be understood from states near Ep, so 
we begin by analyzing the ground state band and magnetic structure of the low-temperature 
(striped) phase of CaFe2As2 within the LSDA. Results for tetragonal and orthorhombic 
striped structures are very similar, suggesting that the slight difference in measured a and 
b lattice parameters plays a minor role in the description of magnetic properties. Experi- 
mental lattice parameters were employed jol. CaFe2As2 has an internal parameter z which 
determines the Fe-As distance R^e-As- Here we treat 2; as a parameter whos main effect is 
to control the magnetic moment M, which in turn strongly affects the character of magnetic 
interaction, as we will show. 

In Fig. [U^a), the density of states N{E) is shown over a 2 eV energy window for a small 
and large moment case (M=0.4/iB and l.l/is). Particularly of note is a sharp peak near 
Ep, of width ~50 meV. This narrow band is found to consist almost entirely of majority- 
spin Fe dxy and dyz orbitals. The peak falls slightly below Ep for M=OAfiB and slightly 
above for M=l.lfiB- Thus, as -Rpe-As is smoothly varied so that M changes continuously, 
this peak passes through Ep. As a result N{Ep) reaching a maximum around M=0.8/is 
(Fig. m^)). This point also coincides with the LSDA minimum-energy value of i?Fe-As- 
That this unusual dependence originates in the magnetic part of the Hamiltonian can be 
verified by repeating the calculation in the nonmagnetic case. As Fig. Wib) shows, the 
nonmagnetic N{Ep) is large, and approximately independent of Rpe-As- In summary, the 
magnetic splitting produces a pseudogap in N{E) for large M; the pseudogap shrinks as 
M decreases and causes a narrow band of Fe d states to pass through Ep, creating a sharp 
maximum in N{Ep) near M=0.8fiB- As we will show below that moment demarcates a 
point of transition from itinerant to localized behavior (see also Ref. jl3l]). 

Next we turn to magnetic excitations. In the standard picture of magnons in metals, 
Imxo(i^) is significant only for frequencies exceeding the magnetic (Stoner) splitting of the d 
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bands, e^^—elj=IM [cf Eq. [2]). In such cases Imxo is small at low q and low energies and well 
defined magnons appear at energies near |1 — /Rexol = (c/ Eq. [1]). As Imxo increases, it 
initially broadens the (formerly sharp) SW spectrum (X'(q); as it becomes large the spectrum 
can become incoherent, or (Stoner) peaks can arise from Im^o, possibly enhanced by small 
1 — /Rexo- Fig. [D^c) shows Imxo(Q'=0, uj) on a broad energy scale. The picture we developed 
for N{E) leads naturally to a classification of Stoner transitions into three main types. 

(1) Excitations between the (largely) Fe d^^ states centered near —0.5 eV to d states 
centered near 1 eV (for M=l.lfiB)- The magnetic splitting of these states matches well 
with the usual Stoner splitting IM in localized magnets, and scales with M. {1^1 eV in 
3d transition metals.) These high-energy transitions are depicted by a large red arrow in 
Fig. lH^a) for M=l.lfiB, and by vertical arrows in Fig. Wic) showing xo(9=0,w) for M=0.4, 
0.8, and Ll/i^. They are too high in energy to be observed by neutron measurements. 

(2) Excitations from d^z to d^y and dyz states just above the psueudogap, depicted by 
a small red arrow in Fig. W^a), and slanting arrows in Fig. Wl^d). These are the excita- 
tions probably detected in neutron measurements [l, 2|. This pseudogap is well defined for 
M>l.lfiB, but is modified as M decreases, which leads to the following: 

(3) Near M=0.8fiB, the narrow d band passes through Ep, opening up channels, not 
previously considered, for low-energy, particle-hole transitions within this band. When M 
reaches 0.4/1^, this band has mostly passed through Ep and the pseudogap practically 
disappears. 

What makes the pnictide systems so unusual is that Imxo is already large at very low 
energies (~10 meV) once the sharp peak in N{E) approaches Ep. One of our central find- 
ings is that this system undergoes a transition from localized to a coexistence of localized 
and itinerant behavior as M decreases from M>l.lfiB to M^O.S^b- Moreover, the itiner- 
ant character is of an unusual type: elementary excitations are mostly single particle-hole 
like: they can be well defined in energy and q. Those represent coherent excitations. The 
dependence of N{Ep) on M is not only responsible for them, but also may explain the 
unusual linear temperature-dependence of paramagnetic susceptibility, and the appearance 
of a Lifshitz transition with Co doping jlo| . 

Fig. [2] focuses on the AFM line, q=[if00]27r/a. Panel (a) shows the full Im^ (w) for 
M=l.lfiB for several g-points spanning the entire line, 0<H<1. At low q, peaks Co in x 
are sharp; and u depends on H in the expected manner (uocH). u reaches a maximum 
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FIG. 2: Imx(q, w) along the AFM axis, q=[ff,0,0]27r/a. (a) lm.x{io) for various and M=l.lfiB- 
Extracted from lmx{H,uj) are peak positions (jj{H) (meV) and HWHM T. (b) and (d) depict u 
and the ratio r/ti), for the three moments shown in the key. Panel (b) also depicts, as solid lines, 
contours |1 — 7Rexo| = in the (w, H) plane, for M=l.l and ^A^b- Panel (c) shows lva.x{H=Q.9,uj) 
for the three moments. 

near H=l/2 (Fig. [21^6)), for all three values of M. The peaks broaden with increasing 
g; nevertheless we can associate them with magnons, because they coincide closely with 
vanishing |1 — JRexol and F is not too large. The magnon character is preserved for most 
q at all moments, as Fig. Wyb) shows: but for H>0.8 the peaks are strongly broadened, 
especially for small M. u is in good agreement with neutron data of Ref.|2], except that 
neutron data are apparently smaller than large M calculations predict 

One experimental measure of the validity of the local-moment picture is the ratio of half- 
width at half- maximum (HWHM) F to u, shown in Fig. ^d). For large M and most of q, 
r/u ~0. 15-0. 18. For intermediate and small M, T/u >0.2 for a wide diapason of q. This 
is significant: it reflects the increasing Stoner character of the elementary excitations. Were 
there an abrupt transition into a conventional Stoner continuum as argued in Ref. it 
would be marked by an abrupt change in T/uj. This is not observed; yet damping appears 
to increase with energy and q, reflecting normal metallic behavior. 

Spectra for H=0.9 (Fig. [2]^c)) adumbrate two important findings of this work. When 
M=0.8fiB, a sharp peak in xi^) appears near 10 meV. There is a sharp peak in Xoi^) at 
uj^lO meV also, classifying this as a particle-hole excitation originating from the narrow d 
band depicted in Fig. [H Being well defined in energy it is coherent, analogous to a SW, 
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only with a much larger T /uj. Yet it is strongly enhanced by collective interactions, since 
|1 — /Rexol ranges between 0.1 and 0.2 for a;<100 meV. This new kind of itinerant excitation 
will be seen at many values g, typically at small q. The reader many note the sharp rise 
and fall in T /u at small q in Fig. Wid)- This anomaly reflects a point where a SW and a 
particle-hole excitation coalesce to the same u. 

Returning to -?/=0.9 when M=l.l, there is a standard (broadened) SW at 200 meV; cf 
contours in Fig. Wyb). A second, low-energy excitation can be resolved near 120 meV. This 
is no corresponding zero in the denominator in Eq. ([1]), but no strong peak in Xf){'^)) either. 
This excitation must be classified as a hybrid intermediate between Stoner excitations and 
SWs. Only a single peak remains when M=0.4/xs and the peak in N{E) has mostly passed 
through Ep (Fig. [T]). 
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FIG. 3: Imx(q,w) along the FM axis, q=[0,K,0]27r/a. (a) Imx(u;) at K= 1/12, 3/12, 5/12, 7/12, 
9/12, 11/12, for M=0.8fXB- Strong Stoner excitations can be seen for i^=l/12 and 3/12 near 
cD~10-20 meV. (b): Contours |1 — IRexol = in the {u},K) plane, for AI=1.1ij,b and O.SfiB, 
analogous to Fig. and dominant peak positions w obtained by a nonlinear least-squares fit of 
one or two gaussian functions to xi^) over the region where peaks occur. 



Along the FM line q^={0,K,0)2'n'/a, x{^) is more complex, and more difficult to interpret. 
At M=l.lfiB sharp, well defined collective excitations are found at low q, and broaden with 
increasing q. There is a reasonably close correspondence with the zeros of |1 — JRexol and 
peaks in x, as Fig.[3](&) shows. The M=0.8fiB case is roughly similar, except that for K>0.4: 
excitations cannot be described by a single peak. Note that for fixed q, propagating spin 
fluctuations (characterized by peaks at u) can exist at multiple energies — the magnetic 
analog of the dielectric function passing through zero and sustaining plasmons at multiple 
energies. Peaks in x{^) can broaden as a consequence of this; the K=7/12, 9/12 and 11/12 
data of Fig. [3](a) are broadened in part by this mechanism as distinct from the usual one 
(intermixing of Stoner excitations). Second, consider how du/dq changes with M for K>0.5 
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(Fig. ^b)). When M=l.lfiB, ^ increases monotonically with K. For M=0.8fiB, ^ has a 
complex structure but apparently reaches a maximum before K reaches 1. That du/dq 
changes sign is significant: it marks the disappearance of magnetic frustration between the 
ferromagnetically aligned spins at low moments, and the emergence of stable FM order along 
[010], characteristic of local-moment behavior (see also Ref. 12[). Experimentally, Ref. 2| 
reports du}/dq>0 for K>0.5. While our calculations provide a clear physical interpretation, 
we note significant differences in the moment where we observe this effect (M=0.8-l.l/iB) 
and the effective spin (3=0.2) used in Ref.j2|. 



(b) 
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FIG. 4: Imx(q,w) along the c axis, q=[0,0,L]27r/c. (a) Imx(w) at values of L listed in the key, 
for M=0.8fiB- When M=1.1ij,b (not shown) the spectrum is well characterized by a single sharp 
peak at any L. When M=0.8 or OAfiB the SW peak remains, but a low-energy Stoner excitation 
coexists with it. The latter are most pronounced for larger L but can be resolved at every L. Peak 
positions (single for M=l.lfj,B, double for M=0.8 and OAfiB) are shown in (b). 



Elementary excitations along the c axis, q=(0,0,L)27r/c, bring into highest relief the 
transition from pure local-moment behavior (collective excitations), to one where coherent 
itinerant Stoner and collective excitations coexist. Collective excitations are found for all 
L and all M. For M=l.lfiB, a single peak is found: excitations are well described by the 
Heisenberg model with weak damping. Comparing Fig. Wl^b) to Figs. W^b) and Wl^b), it 
is apparent that u rises much more slowly along L than along H or K, confirming that 
interplane interactions are weak. When M drops to O.SfiB, a second-low energy peak at 
u' emerges at energies below 20 meV, for small q along [0,7^^,0], and for ah q along [0,0,L], 
coexisting with the collective excitation. Why these transitions are absent for M=l.lfiB and 
are so strong at M=0.8fiB can be understood in terms of roughly cylindrical Fermi surface 
at A;=(1/2,0,A;^). Single spin-fiip transitions between occupied states ej^ and an unoccupied 
states e]^_^q separated by ~10 meV are responsible for this peak (see Eq. |2]). They originate 
from "hot spots" where N^ejj and iV(ek+q) both large. 
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From the SW velocities, {duj / dq)q=o, anisotropy of the exchange couphngs can be de- 
termined We find strong in-plane and out-of-plane anisotropics, and pre- 
dict Vb/va=0.55 and find Vc/va=0.35, where fa=490meV-A. Neutron scattering experi- 



ments 
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151 have measured with Vc/vn^b to be ~0.2-0.5. 

n n 

In summary, we broadly confirm the experimental findings of Refs. [1|, |2[, that there is 
a spectrum of magnetic excitations of the striped phase of CaFe2As2, which for the most 
part are weakly damped at small q and more strongly damped at large q. A new kind of 
excitation was found, which originates from single particle-hole transitions within a narrow 
band of states near Ep, renormalized by a small denominator |1 — /Rexol (Eq. ([I])). They 
appear when the Fe moment falls below a threshold, at which point the narrow band passes 
through Ep. The character of itineracy is novel: excitations occur at low-energy, at energy 
scales typically below T/v, and can be sharply peaked, which should make them accessible 
to experimental studies. The distinction between the two kinds of excitations is particularly 
observable in the anomalous dependence of F/a) on g in Fig. M^d). Collective spin- wave-like 
excitations also are unusual: multiple branches are found at some q. Finally, at higher q 
we find that the SW velocity changes sign, which is also observed experimentally at similar 
values of M. 

Overall a picture emerges where localized and itinerant magnetic carriers coexist and 
infiuence each other. This description falls well outside the framework of a local-moments 
model such as the Heisenberg hamiltonian. The low-energy particle-hole excitations are 
strongly affected by Fe-As distance (and thus by lattice vibrations) while the SWs are much 
less so. While phonons and Stoner excitations are generally considered separately, these 
findings suggest that they strongly infiuence each other. 

This work was supported by ONR, grant N00014-07- 1-0479 and by DOE contract DE- 
FG02-06ER46302. Work at the Ames Laboratory was supported by DOE Basic Energy 
Sciences, Contract No. DE-AC02-07CH11358. 
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